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Studies on the unfolding of the atmospheric neutrino spectrum with IceCube 59 using the 
TRUEE algorithm 

The IceCube collaboration 1 
1 See special section in these proceedings 



Abstract: The measurement of the atmospheric neutrino energy spectrum provides information about the diffuse 
neutrino flux from extragalactic sources. A relative increase of the spectrum toward higher energies could be evidence 
for neutrino producing hadronic processes in the cosmic high energy accelerators, such as active galactic nuclei or 
gamma ray bursts. IceCube is currently the largest neutrino detector on Earth and is placed in the antarctic ice at the 
geographic South Pole. IceCube permits the detection of neutrinos with energies beyond 10 6 GeV. Since the acceptance 
and the resolution of neutrino telescopes suffer from the finite resolution and limited acceptance, a regularized unfolding 
method is used to extract the energy distribution of neutrinos from the measured observables. For AMANDA, the 
unfolding was done with the THAN algorithm. Based on the basic concept of this program and for data analyses in the 
ROOT frame, a new deconvolution algorithm (TRUEE) has been written and tested. With this new algorithm, stud- 
ies on the analysis of the atmospheric neutrino spectrum measured with the IceCube 59 string configuration are presented. 
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Ruhe 2 (tim.ruhe@udo.edu) 

2 Department of physics, TU Dortmund University, D-44221 Dortmund, Germany 

Keywords: IceCube; TRUEE; RUN; regularized unfolding; atmospheric neutrino; energy spectrum 



1 Introduction 

IceCube is the largest neutrino detector ever built and is 
located at the geographic South Pole. It consists of 5160 
digital optical modules (DOM) arranged along 86 strings 
forming a three-dimensional grid covering a cubic kilome- 
ter in the glacial ice [1]. While traveling through the ice the 
high energy neutrino-induced muons produce Cherenkov 
light which can be detected by the DOMs providing direc- 
tional and energy information of the muon track. One of 
the main goals of IceCube is the detection of extragalac- 
tic neutrinos for understanding of cosmic ray production in 
cosmic accelerators. Neutrinos from interactions of cosmic 
rays with the Earth's atmosphere represent a background 
for the extragalactic neutrinos. Thus, a precise measure- 
ment of the atmospheric neutrino flux is important for un- 
derstanding this background. Since the spectral index of 
the flux distribution depending on neutrino energy is lower 
for extragalactic neutrinos (following the spectral behav- 
ior of Fermi accelerated cosmic rays 7^2 [2]) than for 
atmospheric neutrinos (7 <~ 3.7) [3], a contribution of ex- 
tragalactic neutrinos would cause an enhancement of the 
flux in the high energy region of the spectrum. 

The energy of the primary particles is convoluted with the 
interaction probability and the detector finite resolution and 
limited acceptance. Therefore the neutrino energy has to be 



estimated from energy-correlated, measured observables. 
For this purpose a regularized unfolding algorithm is de- 
veloped and applied. 

In this paper the atmospheric neutrino sample from the 
measurement with the IceCube 59 (IC 59) string configu- 
ration is used. The energy spectrum is unfolded with the 
new deconvolution algorithm TRUEE. 

2 Regularized unfolding 

The convolution of the neutrino energy with the interaction 
probability and detector response gives us the measured ob- 
servables in the detector and this relation can be expressed 
as a Fredholm integral equation of the first kind if neglect- 
ing background. From discretization a linear matrix equa- 
tion can be obtained where the measured distribution is a 
product of the detector response matrix and the neutrino 
energy distribution. The response matrix is obtained from 
Monte Carlo (MC) simulation. Thus, an unfolding algo- 
rithm needs as input MC-simulated assumed energy distri- 
bution with the resulting distributions of measured observ- 
ables to determine the detector response and the measured 
observables distributions from data to estimate the neutrino 
energy flux. 
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Every unfolding requires an a-priori assumption about 
some properties of the result realized in a regularization to 
reduce strong negative correlations between unfolded data 
points. Here the Tikhonov regularization [5] is used achiev- 
ing a smooth distribution by minimizing the curvature of 
the result during the unfolding fit. 

3 Unfolding algorithm TRUEE 



4.1 Selection of observables 

As a first step the selection of energy-dependent observ- 
ables is made. TRUEE automatically produces scatter plots 
of the sought-after variable and observables and the re- 
lated profile histograms to check if a correlation is present. 
The inspection of scatter plots with different observables 
showed the correlation between primary neutrino energy 
and the following observables: 



In AMANDA, the precursor of IceCube, the Regularized 
UNfolding (IZUAf) algorithm [6] was used to unfold lepton 
energy spectra [7]. IZUAf was developed in 1984 using the 
programming language FORTRAN 77. Therefore IZUAf is 
not easy to install and use in combination with modern soft- 
ware. TRUEE - Time-dependent Regularized Unfolding 
for Economics and Engineering problems or just TRUE 
Energy is a new software including the IZUAf -based un- 
folding algorithm. 

One special property of the IZUAf algorithm is the 
parametrization of unfolded distribution using a superpo- 
sition of cubic basis splines. The spline coefficients are de- 
termined from the unfolding and the superposed function is 
transformed to the final histogram. At the spline overlaping 
points (knots) the function is continuously differentiable up 
to the second derivative, so that a Tikhonov regularization 
with the second derivatives in the smoothing operator can 
be performed. The user determines the number of splines 
by defining the number of knots. Regularization is con- 
trolled by the number of degrees of freedom, which can be 
given by the user but can also be suggested by the software. 
A small number leads to strong regularization. 

TRUEE has been developed within the Collaborative Re- 
search Centre SFB 823 in Dortmund. Besides the core reg- 
ularized unfolding from IZUAf the new software contains 
user friendly functions, which make the procedure of an 
unfolding analysis more comfortable. The functions used 
in this analysis are explained in the next section using the 
IC 59 neutrino sample. 

4 Unfolding analysis 

For this analysis 10 % of the measured data were used. Af- 
ter event selection a sample of 3160 neutrino events mea- 
sured within the zenith angle range of 88° to 180° with 
IC 59 is obtained. Thus, most of the events in the sample 
were caused by neutrinos having traveled through the Earth 
before undergoing an interaction inside or in the vicinity 
of IceCube. For 100% of the data we expect more than 
30000 neutrino-induced events and thus higher statistics in 
the high energy region. Based on a Monte Carlo study the 
purity of the sample is estimated to be higher than 95 %, 
therefore the background formed by atmospheric muons 
is neglected. The sample is obtained using straight pre- 
cuts followed by an event selection using the multivariate 
method Random Forest [8] within the framework Rapid 
Miner [9]. 



• Number of DOMs having a signal 

• Number of strings with at least one hit-DOM 

• Track length in a certain time window (MPE- 
FitJJDirC) 

Figure 1 shows the correlation between neutrino energy and 
number of DOMs. 




1 1.2 1.4 1.1 



2 2.2 2.4 2.6 

log(number of DOMs) 




2 2.2 2.4 2.6 

logfnumber of DOMs) 



Figure 1 : An example of scatter plot and related profile his- 
togram to check the correlation between the energy and the 
observable (here number of DOMs). An optimal correla- 
tion is present in a monotonically changing profile function 
with small uncertainties. 

Since TRUEE is able to use up to three observables at the 
same time for unfolding, different binnings of the selected 
observables for the response matrix have been checked. 
The most suitable binnings could be chosen by inspecting 
results after running unfolding in the test mode (see Sec. 
4.2). 
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4.2 Test mode 

An unfolding algorithm expects input from the user con- 
cerning some parameters such as the number of bins for 
histograms or the degree of regularization. To check which 
user-defined parameter settings give the optimal unfolding 
result a test mode is included in TRUEE. In this mode only 
simulations are used. Since we neglect the atmospheric 
muon background, the MC sample contains only neutrino 
events after application of event selection techniques. The 
energy distribution of simulated neutrino events has been 
reweighted, so that the generated flux follows the atmo- 
spheric neutrino flux predicted by Honda [3] containing 
a prompt component from Naumov Recombination Quark 
Parton Model (Naumov RQPM) [4]. The prompt com- 
ponent consists of neutrinos from decays of short-lived 
mesons containing charm quarks. 

An MC sample that is statistically equivalent to the ex- 
pected experimental data sample is used as a pseudo real 
data sample for unfolding. Since the real sought-after dis- 
tribution is known in this case, it can be compared to the 
unfolded distribution. The optimal parameter settings are 
chosen with an L-curve approach [10] by examining the 
tradeoff between regularization strength and fit to the true 
distribution. For the IC59 sample the following parameter 
settings are used for the final unfolding: 

• Number of DOMs with 20 bins 

• Number of strings with 10 bins 

• Track length with 10 bins 

• Number of knots 16 

• Number of degrees of freedom 5 

The unfolded test spectrum is shown in Fig. 2. 
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Figure 2: Unfolding of pseudo data compared to the real 
distribution using test mode. The unfolding does not con- 
sider the interaction probability, acceptance or systematic 
uncertainties. 



4.3 Unfolding result 

The unfolding procedure with the parameter settings deter- 
mined in Sec. 4.2 can now be applied to the IC 59 neutrino 
sample. The generated MC neutrino sample for determina- 
tion of the detector response contains only simulated events 
that undergo an interaction within or close to the detector. 
This procedure is necessary to reduce simulation time and 
memory. After passing all event selection steps the final 
sample contains only a fraction of neutrino events. Thus, 
the unfolded distribution represents only neutrinos which 
interacted, triggered the detector and passed the event se- 
lection (Fig. 3). 
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Figure 3: Unfolding of IC59 data gives the distribution 
of selected neutrino events depending on energy. The un- 
folding does not consider the interaction probability, accep- 
tance or systematic uncertainties, yet. Furthermore with the 
full data sample the number of events per bin will increase 
by a factor of ten. 

To calculate the neutrino flux for all neutrinos within the 
zenith angle range, the unfolded spectrum has to be di- 
vided by the effective area. The effective area is the ratio 
of observed event rate and incoming flux and depends on 
the properties of the selected event sample. It includes the 
muon neutrino cross section, the probability for the muon 
to be detected and the detector efficiency for muon detec- 
tion and event reconstruction. The effective area for the 
current sample is shown in Fig. 4. 

The effective area is rising at higher energies due to the in- 
creasing cross section of neutrinos and to the higher length 
of the muon tracks. Therefore the probability to detect and 
reconstruct such a long track is rising. For the events with 
vertical upgoing tracks the effective area is decreasing be- 
cause of the rising probability for absorption of neutrinos 
by the Earth. 

We demonstrate the performance of the unfolding tech- 
nique by showing an example in Fig. 5 of how an un- 
folded energy distribution (Fig. 3) can be translated into 
a neutrino flux spectrum when the effective area (Fig. 4) 
is known. Additional spectra have been unfolded with the 
same parameter settings but with different assumptions of 
the neutrino flux in MC to check the possible bias intro- 
duced by the assumption. Shown are results trained with 
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Figure 4: Effective area for the current neutrino sample de- 
pendent on neutrino energy. Illustrated are areas for dif- 
ferent zenith angle ranges and for the average of the whole 
zenith range considered in the analysis. 

MC weighted to atmospheric (Honda), to atmospheric with 
prompt (Honda-Naumov) and to atmospheric with prompt 
and 1.6 • 10~ 8 -E~ 2 neutrino flux. All three results have 
only small deviation in the low statistics region, thus the 
introduced bias is negligible. 
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Figure 5: An example of the atmospheric neutrino energy 
spectrum from 10 % of IC 59 data unfolded with TRUEE. 
Shown are three unfolding results, using different MC dis- 
tributions to determine the detector response. The uncer- 
tainties take into account statistics and bin-to-bin correla- 
tion, determined by the unfolding software. The spectrum 
is weighted by squared energy for a better illustration. 

Since a continuous function is unfolded taking into account 
event migration between bins, the last bin is estimated to be 
non-zero even though the statistics of current sample is low. 

4.4 Verification of simulation 

A function to verify the result was developed in VMM and 
was transferred to TRUEE. The user has the possibility to 
check whether the simulation of all observables agrees with 
the experimental data and thus verify the unfolding result. 
After the unfolding, the MC events are reweighted by the 
unfolded distribution. The MC sample describes the real 



data now. In this case all observables, not only those which 
have been used for the unfolding, should match in their 
distribution the measured data. The observables whose 
distributions do not match are not correctly simulated. If 
none of the distributions match, probably the unfolding did 
not work properly. In this case the user should check if 
the simulation of the detector response was correct. For 
the unfolding example shown in this paper the verification 
showed an agreement between the experimental data and 
reweighted simulations. 

5 Conclusion and outlook 

The new unfolding algorithm TRUEE shows a good perfor- 
mance in estimation of an atmospheric neutrino spectrum. 
The algorithm is able to estimate a steep distribution cover- 
ing several orders of magnitude and thus is a dedicated tool 
for astroparticle physics. The analysis is facilitated by ad- 
ditional functions and an easy installation and ease of use 
of the software. 

A new energy region of the atmospheric neutrino flux can 
be explored with IC 59. The simulation predicts an exten- 
sion to energies up to 10 6 GeV. The estimation of the neu- 
trino energy spectrum with 10% of the IC59 data will be 
done by unfolding with TRUEE. The subsequent unfolding 
of the 100 % IC 59 data sample is expected to determine the 
energy spectrum with more precision in the high energy re- 
gion due to higher statistics. 
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Abstract: One of the guaranteed fluxes under study by the IceCube neutrino telescope are neutrinos originating from 
cosmic ray induced air showers. These neutrinos come from the decay of n and K mesons (the conventional flux) and 
from the decay of charmed mesons (the prompt flux). Although several flux predictions exist, the electron neutrino flux 
has been measured only up to GeV energies. At TeV energies, where atmospheric neutrinos are an inevitable source 
of background events for astrophysical neutrino searches, the prompt flux becomes important and the flux predictions 
vary greatly. The detection of electromagnetic and hadronic particle showers, which are not only produced by electron 
neutrinos but which can be found in the final states of charged and neutral current interactions of all neutrino flavours, 
remains challenging. Given the sensitivity to all neutrino flavours, the good energy resolution that will be possible with 
fully contained shower events and the possibility to isolate the prompt from the conventional flux, the prospects of this 
detection channel are very promising. This poster will present an analysis done on a data sample collected with IceCube 
in its 40 string configuration as it was running from 2008 to 2009. The development of the event selection on a small part 
of the sample will be discussed. 
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1 Observing Neutrinos at the South Pole 

The possibility to measure or constrain the flux of astro- 
physical neutrinos could help to solve a number of ques- 
tions of which one of the most prominent, the question 
of the origin of cosmic rays, remains unanswered nearly 
a century after their discovery. Experiments that aim at the 
detection of these neutrinos must compensate for the small 
interaction cross sections and the low expected fluxes with 
increased size. With this year's completion of IceCube [1], 
the biggest neutrino detector to date, such an experiment 
is now available. For the experiment a cubic kilometer of 
glacial ice below the geographical South Pole was instru- 
mented with photomultiplier tubes in order to detect the 
Cherenkov light of charged secondaries generated in neu- 
trino interactions. 

During the last 7 austral summers 86 holes were melted 
2.5 km deep into the ice and into each a cable holding 60 
so-called Digital Optical Modules (DOMs) has been de- 
ployed. The light sensors on 78 of these strings form a 
grid with a horizontal spacing of 125 m and a vertical dis- 
tance of 17 m. As the spacing basically determines the en- 
ergy threshold, the detector center was augmented with the 
denser DeepCore infill array between 2009 and 2010. The 
data-taking started already during the construction phase. 



This work uses data recorded between April 2008 and May 
2009 when 40 strings were operational (IceCube 40). 

IceCube's main physics goal is the detection of astrophys- 
ical neutrinos at energies above 100 GeV. These neutrinos 
must be isolated from the much larger flux of leptons cre- 
ated in cosmic ray induced air showers [2], Among these 
a huge number of muons originating mostly from pion and 
kaon decays form the biggest part of the background. In 
the same air showers also atmospheric neutrinos are cre- 
ated [3]. In order to separate them from the astrophysical 
neutrinos a good understanding of their energy spectrum, 
flavour ratios and angular distribution will be helpful. This 
in turn is tightly coupled to our knowledge of the cosmic 
ray composition and hadronic interactions at energies that 
are out of reach of accelerator experiments. 

The atmospheric neutrino spectrum is expected to consist 
of two components, the conventional flux from decaying 
pions and kaons [4, 5] and the prompt neutrinos from de- 
cays of short lived charmed mesons [6, 7]. The existing 
flux predictions for the latter vary widely and current mea- 
surements of the muon neutrino flux [8] are not yet able 
to resolve any prompt from the conventional component 
(see Fig. 1). Compared to muon neutrinos the flux of at- 
mospheric electron neutrinos is lower and falls with a sim- 
ilar steep power law. Taking advantage of the lower en- 
ergy threshold of the DeepCore array, IceCube has recently 
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Figure 1 : Different flux predictions for atmospheric neutri- 
nos (taken from [4]-[7]) compared to a measurement of the 
atmospheric i^-flux with IceCube 40 [8]. 

detected atmospheric neutrino induced showers around a 
mean energy of 40GeV [9]. However, at TeV energies 
this measurement remains challenging, and only recently 
an analysis on the same IceCube 40 dataset started to find 
several promising candidate events [10]. In the energy 
spectrum of neutrino induced particle showers the prompt 
component is expected to emerge from the conventional at 
about 10 5 GeV which is about an order of magnitude lower 
than for muon neutrinos (see Fig. 1 and [11]). This makes 
shower events a suitable tool to isolate the prompt compo- 
nent. 

2 Neutrino Induced Particle Showers 

The events of interest in this study are particle showers 
emerging from deep-inelastic neutrino nucleon scattering. 
Particle showers can be found in the final states of charged 
current (CC) electron and tau neutrino interactions and in 
all neutral current (NC) interactions. Since IceCube can- 
not distinguish u e and low-energy v T CC interactions from 
all-flavour NC interactions, analyses tailored to this event 
signature are effectively sensitive to all neutrino flavours. 
In NC interactions neutrinos deposit only parts of their en- 
ergy so they show up as less energetic cascades. This leads 
to a lower effective area for muon neutrinos. 

At TeV energies the particle showers have lengths of a few 
meters. But due to the large DOM spacing and the scatter- 
ing of light in the ice showers appear as nearly point-like 
light sources. This results in spherical hit patterns which at 
higher energies appear significantly different from the hit 
patterns of muon tracks. 

The separation from the muonic background is mostly im- 
peded by the fact that high energetic muons stochasti- 
cally undergo catastrophic energy losses in the form of 



bremsstrahlung showers along the track. Because of the 
considerable energy deposition these bright electromag- 
netic showers change the appearance of the track and make 
them less distinguishable from the searched signal. This 
has also a connection to the cosmic ray composition be- 
cause proton air showers produce more often single highly 
energetic muons than for e.g. iron showers. From the lat- 
ter often whole bundles of muons reach the detector and 
traverse the detector nearly in parallel. As the individual 
muons will have their stochastic energy losses at different 
positions, the whole bundle appears sufficiently different 
from a single particle shower and is easier to reject. Exten- 
sive simulations performed in the context of similar anal- 
yses done on the IceCube 22 dataset confirmed this effect 
albeit with low statistics. Those muons which passed all 
cuts were originating from proton air showers [12]. 

For electromagnetic showers the light yield scales lin- 
early with energy. It has been shown in a Monte Carlo 
study that for electron neutrino interactions with energies 
of lOTeV-lPeV and well contained interaction vertices 
the energy may be reconstructed with a precision of about 
Alog 10 (^) = 0.13[17]. 

3 Event Selection 

In order to minimize statistical bias a blind analysis is per- 
formed. From the 364 days of usable data, 32 days are 
chosen to develop the event selection. The data was sam- 
pled uniformly over the year in order to reflect seasonal 
variations in the muonic background rate. Secondly, a 
large background sample of simulated muons from more 
than 10 12 air showers were generated. A version of COR- 
SIKA [13] with the Sibyll interaction model was adapted 
for IceCube and used to simulate the Horandel polygonato 
cosmic ray spectrum [14]. Additionally more statistics of 
protons are currently produced in order to study the impact 
of composition uncertainties on the background prediction. 
For the expected signal interactions electron, muon and tau 
neutrinos were generated with a collaboration-internal sim- 
ulation package that is based on ANIS [15]. 

The IceCube 40 detector operated at a trigger rate of about 
1300 Hz. An online event selection based on two quickly 
calculable variables selected events at a rate of 16 Hz. A 
straight line fit through all hit DOMs at position x^ and 
time ti yields a parametrization of the form Xi = x + 
v ■ (U — to) where the parameter \v\ denotes how fast the 
hit pattern evolves. The second variable uses an analogy 
to classical mechanics in which it interprets the hit pattern 
as a rigid body and the recorded amount of light in each 
DOM as a mass. Spherical hit patterns can then be selected 
by calculating the eigenvalues of the tensor of inertia and 
requiring that all three eigenvalues are nearly of the same 
size. This online filter starts to get efficient above an energy 
threshold of about 1 TeV and is optimized for the search for 
the expected astrophysical E~ 2 flux for which it yields an 
efficiency of about 73%. For the less energetic atmospheric 
electron neutrino flux the efficiency is only about 35%. 



32nd International Cosmic Ray Conference, Beijing 201 1 



The selected events were transmitted via satellite to insti- 
tutes in the North where more elaborate track and vertex 
likelihood reconstructions can be performed. They provide 
a sufficient angular resolution for incident muons and with 
the likelihood value of the vertex reconstruction a quality 
parameter to select particle showers. Based on the vertex 
reconstruction also an energy estimator that considers the 
depth dependent optical properties of the ice [16] is run. 
Cuts based on these variables reduce the data rate to 2 Hz 
while keeping about 60% of the atmospheric v e signal. Ac- 
cording to a predicted atmospheric neutrino flux [5] the 
sample contains at this point about 1200 v e and 10000 
(CC+NC) events which are still buried below 50 • 10 6 at- 
mospheric muon events. The effective areas up to this cut 
level are shown in Fig. 2. 

All passing events are fed into a more elaborate likelihood 
reconstruction [17]. This takes into account the full timing 
and amplitude information of the recorded light as well as 
tabulated results of detailed simulations of how light prop- 
agates in the ice [18]. For showers this provides estimates 
for the time and position of the interaction as well as the 
amount of deposited energy. Also the track reconstruction 
is repeated with an iterative optimization strategy in order 
to avoid local minima of the likelihood and to improve the 
angular resolution for background events [19]. 

So-called split reconstructions, which split the recorded 
photons by time into two sets and reconstruct each set indi- 
vidually, provide further information about the event due to 
the different timing behaviour of tracks and showers. For 
a track, later hits are downstream along the track while for 
particle showers they are centered around the vertex but at 
larger distance. 

Based on an argument that shower induced hit patterns 
should be spherical another cut variable can be constructed. 
For an imaginary sphere with a given radius and centered 
at the reconstructed vertex one can calculate the fill ra- 
tio A J hit i M J sphere, where A 'sphere denotes the number of 
all DOMs in the sphere and Nhu the number of triggered 
DOMs. This is especially useful to reject events containing 
several coincident atmospheric muons because the hit pat- 
terns of e.g. two coincident but well separated tracks can 
have many untriggered DOMs in between. 

In order to further reduce the muonic background, DOMs 
at the surface of the instrumented volume are used to veto 
tracks that appear to enter and traverse the detector. Events 
for which the first triggered light sensor is located on the 
outer layer of the detector are rejected. Together with the 
requirement that the reconstructed vertex is located inside 
the fiducial volume, this forms a strict containment cut. 

Finally those variables that still provide separation power 
are combined with a machine learning algorithm. In the 
TMVA framework [20] a boosted decision tree is trained 
which provides a final cut variable to select particle show- 
ers. Current investigations suggest that this event selection 
is able to remove the remaining background while keep- 
ing the prospect to find a few atmospheric neutrinos in the 
whole sample. 
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Figure 2: Effective areas for electron and muon neutrinos 
for the early stages of the analysis. The width of the bands 
denotes the statistical error. This differs between cut levels 
because datasets of different size have been used. The peak 
at 6.3 PeV for electron neutrinos is due to the Glashow res- 
onance. The drop in effective area for muon neutrinos be- 
tween trigger and online filter level illustrates the effect of 
tailoring the analysis to neutrino induced particle showers. 
At the presented cut levels the contribution of muon tracks 
from charged current interactions is still present. Therefore 
the effective area is still higher for muon than for electron 
neutrinos. 

However, this statement relies on the Monte Carlo back- 
ground prediction which has to be scrutinized before un- 
winding. Accordingly, studies of the systematic uncertain- 
ties in the simulated background sample (like for example 
a lighter cosmic ray composition) are ongoing and will be 
presented together with the final event selection at the con- 
ference. 
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Abstract: For the first time atmospheric neutrino-induced cascades are observed in IceCube, using the 79-string detector 
with the DeepCore extension. Using 23 days of test data, a high statistics sample obtains 824 events with a cascade 
expectation of 312 events, and using 281 days of full data, with tighter cuts 1029 events are observed of which 59 % are 
predicted to be cascades. 
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1 Introduction 

The main mission of the DeepCore extension to IceCube, 
shown in Fig. 1, is to observe neutrinos down to energies 
as low as 10 GeV [1]. At these energies there are valu- 
able science topics such as neutrino oscillations [2] and 
neutrinos from low mass WIMP dark matter annihilation 
[3]. DeepCore uses denser module spacing, upgraded pho- 
tomultiplier tubes (PMTs), and deployment in the clearest 
ice along with a lower trigger threshold than the surround- 
ing IceCube detector [4] to enhance low energy neutrino 
detection. 

In this paper, we discuss the performance of the first year 
DeepCore with 79 strings of IceCube installed (IC-79). 
In this configuration the DeepCore sub-array includes 6 
densely instrumented strings optimized for low energies 
plus the 7 adjacent standard strings. Finally, we present 
preliminary results on the first observation of atmospheric 
neutrino-induced cascades in IceCube. 

2 Performance 

The IC-79 detector finished a year-long data-taking cycle 
in May 201 1 and these data have been processed and re- 
constructed. 

Full PMT waveforms are read out from digital optical mod- 
ules (DOMs) in hard local coincidence (HLC), which re- 
quires both the initial DOM and the immediate neighbor or 
next nearest neighbor DOM to register a hit within a time 
span of ±1000 ns. IceCube also records hits in soft lo- 
cal coincidence (SLC), for which no neighboring hits are 




Figure 1: Schematic top and side views of IceCube in its 
79-string configuration. The DeepCore infill array and the 
relative dust concentration are also shown. 



required [5]. Although the DOM noise rate is a factor of 
50 higher in SLC mode than in HLC mode, with current 
software cleaning algorithms, SLC hits from physics inter- 
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actions can be identified with over 90 % purity. This pro- 
vides an advantage at low energies where the fraction of 
SLC hits is significant, because including the SLC hits im- 
proves reconstruction, background rejection, and particle 
identification. 

A low trigger threshold (SMT3) is applied to all DOMs 
in the fiducial region (shaded area below the dust layer in 
Figure 1) by requiring 3 or more HLC hits within a 2500 ns 
time window. Additionally, the high quantum efficiency 
PMTs, with a 35 % increase in light collection with respect 
to regular IceCube PMTs, help trigger on neutrinos with 
energies as low as 10 GeV. 

3 Observation of Atmospheric Cascades 

Atmospheric neutrinos are the products of the interaction 
or decay of pions and kaons produced when cosmic ray 
primaries interact with nucleons in the upper atmosphere. 
Charged-current (CC) electron and tau neutrino interac- 
tions, and neutral-current (NC) neutrino interactions of any 
flavor, produce cascades that are approximately spherical in 
light emission. However, IceCube can not distinguish be- 
tween a CC cascade and an NC cascade at the considered 
energies. While atmospheric muon neutrinos have been 
observed in large quantities by IceCube [6], atmospheric 
cascades have a substantially lower flux and have not been 
identified in previous IceCube analyses. For current cas- 
cade searches above 1 TeV, there are on-going analyses 
with the 40 string IceCube detector [7, 8]. 

3.1 Background 

The dominant background for the atmospheric cascade 
analysis is penetrating cosmic ray muons. New rejection 
techniques are required to remove muon background in 
large part because the poor directionality in cascades as 
compared to tracks makes it difficult to use the Earth as 
a shield as is done in conventional v^ c detection [6]. The 
signal is therefore required to be contained in DeepCore, 
and veto techniques are applied to remove about six orders 
of magnitude of background events while retaining rea- 
sonable signal efficiency for atmospheric neutrino-induced 
cascades in the fiducial volume. 



3.2 Event Selection 

IC-79 was operational from June 2010 to May 2011 with 
over 90 % up-time for physics analyses. The DeepCore 
SMT3 trigger fired at a rate of 180 Hz, which was reduced 
to ~ 17 Hz using an on-line filter run at the South Pole. In 
the DeepCore on-line filter, an estimate of a neutrino inter- 
action vertex and time is obtained by calculating the center 
of gravity (COG) of all HLC hits in the fiducial region. An 
event is rejected when the velocity between an HLC hit in 
the veto region and the COG is consistent with an incoming 
muon. A 90 % cosmic ray muon background rejection is 
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Figure 2: Output BDT5 Distribution based on the five vari- 
ables. Events with a BDT5 score >0.067 are kept. 



achieved with respect to the trigger while keeping > 99 % 
atmospheric neutrinos that interact in the fiducial volume. 
Noise cleaning algorithms remove hits which are not corre- 
lated in space and time with other hits. Then, events with at 
least 8 total remaining hits and at least 4 hits in the fiducial 
region are selected. 

In order to attain a manageable data volume for time- 
intensive reconstructions, the first background rejection is 
achieved by an initial selection with a Boosted Decision 
Tree (BDT) [9]. The "BDT5" has 5 simple observables 
which quantify the event topology in the detector. These 
are constructed from the position and shape of the hits in 
the detector, and their time and charge. Figure 2 shows the 
output distribution and the reduction cut is made at 0.067. 
After the BDT5, the data rate is reduced to 0.1 Hz and the 
atmospheric v e rate is predicted as 6.3 x 10~ 4 Hz, corre- 
sponding to 63 % retention with respect to the trigger. 

The reduced data set was then processed with iterative 
likelihood reconstructions taking into account detailed 
Cherenkov light propagation in the ice [10]. Hits falling 
outside a time window [ -3000 ns, +2000 ns ] with respect 
to the trigger, or outside of a 150 m radius of a neighboring 
hit within a 750 ns time window are removed. Then, with 
the cleaned hits, after demanding 8 or more hits in the fidu- 
cial region within a 1000 ns sliding window, another BDT 
is formed with 7 input parameters (BDT7). The first two 
variables measure the locations of the earliest hits in terms 
of radial and vertical coordinates to select contained events. 
The next three variables separate cascade-like events from 
muon-like events; an event is split in half then charge de- 
position, COG, and particle speed are compared between 
the two separate halves. The final two compare likelihoods 
of a cascade hypothesis to that of a muon hypothesis. The 
top left plot of Figure 3 shows the output BDT7 distribu- 
tion. A cut at BDT7> 0.22 reduces the atmospheric muon 
background to 5.0 x 10 -4 Hz by rejecting a factor of 200 
more (3.6 x 10 5 cumulatively) background while retaining 
~ 40 % of the v e signal (2.6 x 10~ 4 Hz) compared to the 
previous BDT5 cut. 
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Figure 3: The output distribution of BDT7 in shown upper left corner and extra cuts to select two neutrino samples are 
listed as well. These plots are after BDT5 cut. Solid vertical lines indicate the final Loose Cut sample while the dotted 
vertical lines are for final Hard Cut sample. Arrows indicate events kept. 



3.3 Results 

Two sets of cuts on input variables to BDT7 were used for 
the final selection of neutrino-rich samples. The first aims 
for a high efficiency (Loose Cut) while the second aims for 
a higher purity (Hard Cut). Figure 3 shows two combi- 
nations of these cut options. The Loose Cut option is indi- 
cated in solid vertical lines. Two containment cuts Zi st (up- 
per right) and Z rcco (lower left) based on the vertex depth 
measurements ensure that most signal events are well con- 
tained inside the DeepCore fiducial volume. The NdRdT 
(lower right) extends the on-line veto method by adding 
a travel speed window to catch only muon-produced light 
in the veto region and therefore removes residual muon 
background events. Finally, the cut RLLH casca d c (upper 
middle) selects events that fit a cascade hypothesis well, as 
measured by the log likelihood from a fit. With the Loose 
Cut options, the events observed are consistent with v GC , 
v^ c , and v e interactions with only two cosmic ray muon 
events expected in 28 hours of simulation. 824 events are 
observed in 23 days of test data with 590 total background 
event expectation (v GG + y). Systematic errors are not in- 
cluded. 



Signal simulation of atmospheric neutrinos predicted ap- 
proximately 3800 cascades and 7200 v GG events in 281 
days (full oscillations for and v e are included based on 
the parameters in [11]) as shown in Figure 4 with an esti- 
mated mean energy of 40 GeV. Unfortunately, events below 
10 GeV are not accurately simulated in current neutrino 
generation and are not included in this work. When con- 
sidering oscillations of atmospheric neutrinos the rate of 
v GC events is reduced by 20 % with the Loose Cut selec- 
tion. The effect of oscillations on v e is marginal compared 
to the data rate indicated as a dotted line in Figure 4. The 
v T oscillation from in these energies is not yet included 
and under active development with a new simulation pack- 
age, GENIE [12]. Since there is no distinction among v e , 
v T , and v^ c in the IceCube detector at these energies, a 
small fraction of the 20 % disappeared events may be 
re-introduced as cascades if they oscillate to v T . 

The Hard Cut selection for higher purity is shown in Fig- 
ure 3, with an additional requirement of >19 hit DOMs, 
with the result shown in Figure 5. Compared to the Loose 
Cut option, these cuts are stronger and mainly focused on 
rejecting v GC events. The containment cuts select a vol- 
ume smaller than the nominal DeepCore fiducial region to 
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Figure 4: Loose Cut result. The event rate in Hertz as func- 
tion of months is shown. The squares are 23 days of test 
data while simulations are indicated with the various lines. 
Dotted lines are the rates without oscillations and hori- 
zontal solid lines are with oscillation. The data rate shows 
an excess over the prediction and can be explained with 
two contributions; oscillation and cascades. Here, the sig- 
nal predictions are based on Barr et al. [13] 



Data ( 281 days ) 
Sum of All Monte Carlos 
Atmo. v e +v u NC (Cascades) 
Atmo. v u CC (Tracks) 
Atmo. i (No Events in 28 hours) 




100 120 140 
Number of Hit DOMs 



Figure 5: Hard Cut result. The plot shows an event rate 
as a function of the number of hit DOMs. The sum of all 
Monte Carlo is consistent with 281 days of data rate. The 
cascades are expected to contribute 59 % and the tracks are 
expected to be 41 %. No atmospheric muon background 
events remain in 28 hours. The oscillation effect is less than 
3 % due to higher energies of cascades and not included 
here. 

identify an outgoing track from a interaction. With 
281 days of full data, we observe 1029 total events ( 650 
cascades and 454 tracks are expected from simulation ). 
The remaining v^ c events have short muons with a median 
track length of 80 m where the muon tracks are not easily 
visible due to detector granularity. Half of the cascades are 
predicted to be v e events and the other half are z/^ c events. 
This selection favors higher energy cascades and the mean 
cascade energy increases to 180 GeV, so that oscillation of 
— > v x has a small (< 3 %) effect. The atmospheric 
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824 


312 (590) 






3812 (7208) 




Hard Cut 




N obs 


C sig (Qbg) 


281 days 


1029 


650 (454) 



Table 1: The number of events are shown with two final 
selections. N obs means observed events in real data. C slg 
and C bg refer predictions of the cascade signal and its back- 
ground respectively. Note that the total data (281 days) for 
the Loose Cut sample is not open and the results of the test 
data (23 days) are only shown. 



muon simulation predicts zero events in 28 hours. A sea- 
sonal effect in rates will be studied more carefully with the 
full data. 

3.4 Conclusion 

We observe atmospheric neutrino-induced cascade events 
for the first time with IceCube 79 string detector with the 
DeepCore extension. The preliminary numbers are sum- 
marized in Table 1. Possible sources of systematic er- 
rors include ice modeling, detection efficiency of DOMs, 
neutrino-nucleon cross-sections, and atmospheric neutrino 
flux normalizations. Those uncertainties are under evalua- 
tion. 
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Abstract: The discovery of a cumulative flux of high-energy neutrinos from the sum of all cosmic sources in the Uni- 
verse is one of the central goals of the IceCube experiment. The experimental signature of isotropically distributed 
astrophysical sources is an excess of high-energy neutrinos with a characteristic angular distribution over the background 
of less energetic neutrinos produced when cosmic rays interact with the Earth's atmosphere. Such searches are chal- 
lenging because of systematic uncertainties in these fluxes and the detector response. The distribution of reconstructed 
neutrino energies is analyzed using a likelihood aproach that takes into account these uncertainties and simultaneously 
determines the contribution of an additional diffuse extraterrestrial neutrino component. This analysis is applied to the 
data measured with the IceCube detector in its 40 and 59-string configurations, covering the period from April 2008 to 
May 2010. No evidence for an astrophysical neutrino flux was found in the 40-string analysis. The upper limit obtained 
for the period from April 2008 to May 2009 is d$/dE < 8.9- 10~ 9 GeV" 1 cm" 2 s" 1 sr" 1 at 90% confidence level in 
the energy region between 35TeV and 7PeV. For the 59-string data from May 2009 to May 2010, an improved anal- 
ysis technique including the angular distribution in the likelihood approach is presented. The preliminary sensitivity is 
d$/dE < 7.2-10" 9 GeV- 1 cm" 2 s" 1 sr" 1 . 
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1 Introduction 

The study of cosmic rays is one of the main aspects of cur- 
rent research in astroparticle physics. Despite all efforts, 
charged cosmic rays have not yet revealed their sources. 
A candidate source class is active galactic nuclei, which 
are believed to accelerate particles up to energies of several 
EeV by the mechanism of Fermi acceleration, e.g., in the 
vicinity of their central supermassive black holes. Through 
hadronic interactions with the surrounding matter and ra- 
diation, high-energy neutrinos can be produced. Unlike 
charged cosmic rays and photons, neutrinos propagate al- 
most unaffected by magnetic fields or intervening matter 
through the universe. This makes them an ideal messenger 
particle for astrophysics. 

The neutrino telescope IceCube was built at the geographic 
South Pole with the purpose of detecting neutrinos with 
energies from several tens of GeV to EeV [1]. It consists 
of 86 strings each equipped with 60 optical sensors, dis- 
tributed over an area of roughly 1 km 2 and instrumented in 
depths from 1.5 to 2.5 km in the Antarctic ice. This huge 
volume is necessary to compensate for the very low interac- 
tion probability of neutrinos with matter. After seven years 



of construction, the IceCube telescope was completed in 
December 2010 and is currently the largest detector of its 
kind in the world. 

The detection principle is based on the observation of sec- 
ondary charged leptons and hadrons produced in interac- 
tions of neutrinos in the surrounding ice and rock. These 
emit Cherenkov light which is detected by IceCube's opti- 
cal sensors. From the number of photo-electrons and their 
arrival times, detected by the optical sensors, the neutrino's 
initial direction and energy are reconstructed. Although 
no specific neutrino emitting sources have been discovered 
yet, it is believed that the combined flux of many weak 
sources distributed all over the sky could be detected with 
IceCube. This flux would exceed the flux of cosmic ray in- 
duced atmospheric neutrinos at high energies and would ar- 
rive almost isotropically from all directions. Since it would 
not be possible to identify individual neutrino sources, this 
analysis is known as a search for a diffuse neutrino flux. 

2 Neutrino Event selection 

The first step in the searches for a diffuse astrophysical neu- 
trino flux is to select a sample of neutrino events with high 
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Figure 1: Reconstructed zenith angle distribution of one 
day of experimental data, and of simulated muons from 
air showers and neutrino-induced muons with the 59-string 
configuration at trigger level. The distribution of astrophys- 
ical neutrinos is normalized to the 40-string analysis upper 
limit. 

purity. This contribution presents two searches for a diffuse 
neutrino flux with data from two consecutive years during 
the construction of IceCube. Both analyses focus on the 
selection of high-energy secondary muon tracks. Data was 
taken from April 2008 to May 2009 in the 40-string config- 
uration and from May 2009 to May 2010 with 59 deployed 
strings. The event selections and analysis techniques are 
very similar. The analysis of the 59-string sample has not 
been finalized. 

The reconstructed zenith angle distribution of detected 
events is shown in Fig. 1. The dominant background in 
this analysis are muons from cosmic -ray air showers. At 
trigger level, they outnumber the detected neutrino-induced 
muons by several orders of magnitude. In contrast to neu- 
trinos, muons are easily absorbed by the Earth. Therefore, 
muons produced in the atmosphere enter the detector from 
above and are primarily reconstructed as downward going 
tracks, while muons originating from neutrinos interacting 
with the matter surrounding the detector come from all di- 
rections. 

To reject a large amount of air shower background the 
analysis is restricted to upward reconstructed muon tracks . 
The remaining background is misreconstructed air-shower- 
induced muon tracks, containing a large fraction of muons 
arriving from coincident but independent air showers. For 
the further selection, an algorithm searches for patterns 
separated in space and time in the ensemble of recorded 
light-sensor pulses. This allows rejection of coincident 
events as well as tracks associated with random noise hits. 

For the selection of a high-purity upward-going neutrino 
sample, the remaining data is reduced by a series of quality- 
criteria applied to reconstructed variables like the direc- 
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Figure 2: Distribution of the reconstructed muon energy 
loss for 10% of the 59-string data after neutrino selection. 

tional estimate of the reconstruction. They are described 
in detail in [5]. The final event sample consist of 12877 
neutrino candidate events for the 40-string configuration 
and about 25000 expected events for the 59-string config- 
uration after finalization of the analysis. Based on Monte 
Carlo simulation, the expected contamination of remaining 
background events is less than 1%. 

Figure 2 shows the distribution of the reconstructed average 
energy losses for the selected muon tracks along their path 
in the detector. The experimental data is largely consistent 
with the expectation from atmospheric neutrinos. Most in- 
teresting for this analysis are events with high energy de- 
positions. 

3 Analysis method 

The irreducible background for astrophysical neutrino 
searches consists of conventional atmospheric neutrinos. 
These neutrinos are produced in the decay of pions and 
kaons in cosmic-ray air showers in the Earth's atmosphere. 
They are described by an energy spectrum following a 
power law of about d^/dE oc E~ 3 - 7 and by a character- 
istic zenith angle distribution related to the meson's path 
through the atmosphere. Another - not yet observed - type 
of atmospheric background are so called prompt neutri- 
nos. Prompt neutrinos originate from the decay of heavier 
mesons, typically containing a charm quark[3]. They are 
produced at a higher cosmic -ray energy threshold and be- 
cause of their comparably short lifetimes their energy dis- 
tribution is predicted to follow a harder energy spectrum of 
d^/dE oc E~ 2 - 7 with an almost isotropic angular distribu- 
tion. 

The aim of this analysis is to identify a possible astrophys- 
ical component in the neutrino sample. An astrophysical 
flux can be distinguished from a conventional atmospheric 
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Figure 3: Expected energy (left) and zenith (right) distribution for detected conventional atmospheric neutrinos, prompt 
atmospheric neutrinos and astrophysical neutrinos in the IceCube detector with 59 strings. Left: The astrophysical neu- 
trino flux is normalized to the upper limit (90% CL) of the 40-string analysis presented here. Right: The astrophysical 
and prompt fluxes have been renormalized for better visualization. 



flux by its harder energy spectrum. Assuming shock accel- 
eration in the extragalactic sources, an astrophysical neu- 
trino flux would follow a d<fr/dE oc E~ 2 - power law (see 
Fig. 3). With the presumption of isotropically distributed 
sources over the whole sky, the arrival directions of these 
neutrinos would be isotropic. Their energy spectrum be- 
ing harder than that of conventional atmospheric neutrinos, 
prompt neutrinos are an important background in a search 
for a diffuse flux. 

Relative to the 40-string analysis [5], the ongoing 59-string 
analysis improves the sensitivity to an astrophysical flux by 
considering directional information in addition to energy. 
Figure 3 shows the expected zenith angle distribution of ar- 
rival directions when considering the energy-dependent ab- 
sorption in the Earth, the angular detector acceptance and 
event selection efficiency. The significant differences in an- 
gular distribution for neutrinos of different origin adds sep- 
aration power between the three components. 

A likelihood method is applied to the experimental data to 
fit for the contributions of conventional atmospheric neutri- 
nos, prompt atmospheric neutrinos and astrophysical neu- 
trinos. In the 40-string analysis, the corresponding one- 
dimensional probability density functions (pdf) of the re- 
constructed energy are used to determine the probability 
for an astrophysical and prompt component. For the 59- 
string analysis, two-dimensional pdfs of reconstructed en- 
ergy loss and zenith angle are used to account for both pa- 
rameters and their correlation. Systematic uncertainties are 
taken into account by incorporating nuisance parameters in 
the likelihood function. These uncertainties are discussed 
in the next section. 

The test statistic is a profile likelihood based on a likelihood 
ratio of the best fit of all physics and nuisance parameters to 
the experimental data compared to a fit of only the nuisance 
parameters for each point in the physics parameter space. 



Confidence regions are constructed according to the Feld- 
mann & Cousins approach by generating a large number of 
random experiments based on Monte Carlo simulations [7]. 
In order to estimate the sensitivity of the analysis to a sig- 
nal of diffuse astrophysical neutrinos, random experiments 
assuming the zero-signal hypothesis are generated. 

4 Systematic uncertainties 

A challenge in the search for a diffuse neutrino flux is the 
treatment of systematic uncertainties. Unlike other analy- 
ses of IceCube data, the background cannot be estimated 
from an off-source region in the experimental data. There- 
fore, the background estimation relies on a full-chain de- 
tector simulation. Inputs are, amongst others, air showers 
simulated with CORSIKA [4] and atmospheric neutrinos 
based on [2, 3]. More details can be found in [5, 6]. 

Main uncertainties at high energies are the conventional 
and prompt atmospheric neutrino flux predictions, the cal- 
culated neutrino cross sections and in particular the mod- 
eling of the detector response. Examples for the latter are 
the optical properties of the Antarctic glacial ice and the 
absolute efficiency of the optical sensors. The influence 
of these uncertainties on the final result is determined by 
studying simulations with different settings of these pa- 
rameters. Some uncertainties, such as in the spectral in- 
dex of atmospheric neutrinos, are taken into account with 
nuisance parameters in the likelihood fit. Additional in- 
formation on the sytematic uncertainties can be found in 
[5, 6]. 

One possibly significant uncertainty not taken into account 
in the 40-string analysis is the effect of the knee in the 
cosmic-ray spectrum on the energy spectrum of conven- 
tional atmospheric neutrinos. This leads to an expected 
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Figure 4: Limits and predictions for diffuse muon neutrino fluxes. The thin black lines show the expected flux for at- 
mospheric neutrinos without and with an additional prompt component together with the unfolded atmospheric neutrino 
spectrum by IceCube[6]. The black horizontal lines represent 90%-confidence-level upper limits from different experi- 
ments [8, 9, 5]. The gray curves represent theoretical flux predictions for AGN models [11], gamma ray bursts [12] and 
GZK neutrinos [13]. The thick gray line shows the Waxman-Bahcall upper bound [10]. 



steepening of the neutrino spectrum above several tens of 
Te V, which is within the energy range relevant for this anal- 
ysis but which has so far not been included in our simula- 
tions. The systematic uncertainties related to such a neu- 
trino knee will be incorporated into the 59-string analysis 
using parameterizations of the measured cosmic -ray spec- 
tra. 

5 Results 

The result of the 40-string diffuse neutrino search has been 
published in [5]. The measured energy distribution is con- 
sistent with the expectation from conventional atmospheric 
neutrinos only. No prompt atmospheric nor astrophys- 
ical flux component was found. A small underfluctua- 
tion relative to the expectation was observed in the high- 
energy tail. This results in an upper limit on an astro- 
physical neutrino flux d$/dE oc E~ 2 of d$/dE < 
8.9-10 -9 GeV -1 cm~ 2 s _1 sr _1 integrated over the energy 
range between 35 TeV and 6.9 PeV with 90% confidence. 
This is currently the most constraining limit on a diffuse 
astrophysical neutrino flux and about a factor of two below 
the Waxman-Bahcall upper bound for an astrophysical neu- 
trino flux [10]. At the same confidence level, a prompt at- 
mospheric flux at 70% of the most probable flux predicted 
by Enberg et al. [3] was ruled out. 

The higher statistics of the neutrino sample from the larger 
59-string detector improves the sensitivity to astrophysi- 
cal fluxes by about 35% compared to the 40-string anal- 



ysis. The additional gain from using directional informa- 
tion is about 10% and results in a sensitivity of d^/dE < 
7.2-10" 9 GeV _1 cm" 2 s- 1 sr- 1 (Fig. 4). A further gain of 
roughly 10% in sensitivity for astrophysical fluxes is ex- 
pected when taking into account the effect of the knee in 
the cosmic-ray spectrum. 
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Abstract: IceCube is the first cubic-kilometre scale neutrino observatory dedicated to detecting astrophysical neutrinos. 
A large contribution to the expected neutrino signal is from electromagnetic and hadronic showers (cascades) initiated 
by charged current v e interactions, and neutral current z/ e , Vy. and v T interactions. Cascade energy is reconstructed with 
better resolution than muons, and the atmospheric background is lower. The energy spectrum of astrophysical neutrinos 
is expected to be harder than that of atmospheric neutrinos, so an astrophysical neutrino signal could be observable as a 
break in the cascade energy spectrum. 

Cascades are difficult to detect due to a large background coming from atmospheric muons and muon bundles, many 
orders of magnitude larger than the cascade signal. Large statistics, advanced reconstruction methods and machine 
learning techniques are required to isolate cascade events within the data. There is growing evidence for neutrino-induced 
cascade events in IceCube in several analyses that were carried out using data from April 2008 to May 2009, when 40 
IceCube strings were operational. This is the largest instrumented volume, and yields the most sensitive search for a 
diffuse flux of astrophysical neutrinos using cascades to date. 

Two of the IceCube-40 analyses are described in these proceedings. The high energy analysis measures four background 
events, and sets a 90% confidence level limit for all flavour astrophysical neutrino flux of 9.5 x 10 -8 GeVsr _1 s _1 cm -2 
over the energy range 89 TeV to 21 PeV. The mid energy analysis observes 14 cascade candidate events. The background 
to these events is under investigation. Three of the observed events have reconstructed energies above 100 TeV. 
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1 Introduction 

High energy neutrino production is predicted to occur in re- 
gions of the universe containing astrophysical objects that 
emit large amounts of energy [1]. The same regions are 
predicted to emit the highest energy cosmic rays, whose 
origins are yet unknown. These are dense regions where 
large gravitational forces generate relativistic jets, accel- 
erating charged particles. This is associated with objects 
such as supernovae, gamma ray bursts, and active galactic 
nuclei. High energy neutrinos originating from these ob- 
jects may be observed as a diffuse flux by detectors such as 
IceCube. 

A large proportion of the expected diffuse flux which inter- 
acts with nucleons in the detector results in particle showers 
(cascades). IceCube is capable of detecting cascades pro- 
duced from all flavours of neutrinos. Cascade energy is re- 
constructed with better resolution than that from track-like 
particles such as muons, since cascades are fully contained 
in the detector. Also cascades have lower atmospheric neu- 
trino background flux. The astrophysical flux has a harder 



energy spectrum than that of atmospheric neutrinos, which 
makes diffuse searches a promising route for observing a 
break in the energy spectrum of neutrinos from astrophysi- 
cal sources. A previous cascade analysis in IceCube [2] has 
shown progress towards a detection of atmospheric neutri- 
nos, and set a limit of 3.6x 10~ 7 GeVsr _1 s _1 cm -2 on as- 
trophysical neutrinos (assuming a 1:1:1 flavour ratio) for 
the energy range 24 TeV to 6.6 PeV, using the 22 string 
IceCube detector. There are several cascade analyses us- 
ing the larger 40 string IceCube detector: two low energy 
analyses described in other proceedings at this conference 
[3, 4], and the two high energy analyses described here. 
The goal of these high energy analyses is to search for as- 
trophysical neutrino-induced cascades. 

2 IceCube-40 Data 

IceCube is a Cherenkov neutrino detector located at the 
South Pole. The detector is comprised of Digital Opti- 
cal Modules (DOMs) [5] situated on strings deployed deep 
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in the Antarctic ice. The DOMs house photomultiplier 
tubes (PMTs) [6], along with data acquisition software, 
in a pressure vessel. There are 60 DOMs on each string 
which detect Cherenkov light from charged particles trav- 
eling through the ice. The analyses described here uses 
data from when the detector was still under construction, 
when 40 strings were deployed and acquiring data. Figure 
1 shows the IceCube-40 detector configuration. 
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Figure 1 : IceCube-40 detector configuration. 

The IceCube-40 detector was operational from 6 th April 
2008 to 20 th May 2009. During physics runs all event in- 
formation was sent from DOMs to the surface for process- 
ing if the trigger condition was met. The trigger condition 
requires 8 DOMs to be hit within 5000 ns. The IceCube-40 
physics dataset contains 374 days of data. 

3 Analysis 

The analyses presented here [10, 11] search for an E~ 2 
neutrino flux within the IceCube-40 dataset. These 
searches use cuts on reconstructed event variables, reduc- 
ing the background from atmospheric muons to isolate cas- 
cade events originating from astrophysical sources. 

Neutrino interactions of all flavours were simulated to pre- 
dict the expected signal. This was done using Monte Carlo 
simulations with an energy spectrum of E _1 . These simu- 
lated events were then re-weighted for atmospheric and as- 
trophysical neutrino spectra. The atmospheric re-weighting 
uses the Bartol model [7] for the conventional neutrino flux, 
and the Sarcevic model [8] for the prompt neutrino flux. 
The astrophysical re-weighting uses an E -2 spectrum. 

The dominant background to these analyses comes from at- 
mospheric muons, simulated using Monte Carlo techniques 
with CORSIKA [9]. The simulation was used to train ma- 
chine learning algorithms to develop cuts to separate signal 
from background. 

3.1 Filter Levels of Mid Energy Analysis 

The first level of filtering is run online at the pole in order 
to reduce the data volume to an acceptable level for transfer 



via satellite. The pole filter for the cascade stream during 
IceCube-40 consisted of two cuts placed on reconstructed 
variables. The first variable is a fit to the hit timing, and is 
placed to reject track-like events with a high velocity, and 
keep cascade-like events with low velocities. The second 
variable uses the hit topology of the event defined by Ten- 
sor of Inertia eigenvalues. The cut rejects elongated track- 
like events, and keeps highly spherical cascade-like events 
by cutting on the ratio of lowest eigenvalue to the sum of 
the three. After the data is transferred, level 2 process- 
ing is run, which consists of further reconstructions used 
in higher level filtering. 

The level 3 filter reduces the background further in order 
to run more advanced reconstruction algorithms. This fil- 
ter was applied only to events with a reconstructed cascade 
energy below 10 TeV, where the background is most domi- 
nant. The first cut variable is the reconstructed zenith direc- 
tion assuming a plane-wave track topology. Events coming 
from above the horizon (6* z < 80°) are removed. The sec- 
ond cut variable is the reduced log likelihood. This variable 
is derived from the likelihood that an event is a cascade, 
based on the hit pattern in the detector. Events less likely 
to be cascades are removed. 

The level 4 filter reduces the background further to run 
more sophisticated reconstructions with better variable res- 
olutions. There are three cuts at this level. The first cut 
is on the reconstructed energy, placed at 2.5 TeV. The two 
further cuts at this level are on Spatial Distance and Fill 
Ratio reconstructions. The Spatial Distance cut splits the 
event into two parts based on the timing of hits, with the 
vertex position from each half reconstructed separately. If 
the event is a spherical cascade-like event the two vertex 
positions are expected to be at the same location in the de- 
tector. If the event is an elongated track-like event the two 
vertex positions will be located far apart. This cut requires 
that the two reconstructed vertex positions be within 40 m 
of each other. The Fill Ratio cut defines a sphere based 
on the radius from the mean position of hit DOMs around 
the reconstructed vertex of an event, and looks at the ratio 
of DOMs hit within this sphere over the total number of 
DOMs. If the event is a spherical cascade-like event the fill 
ratio is expected to be close to one. If the event is an elon- 
gated track-like event the fill ratio will be much less than 
one. This cut requires at least 40% of DOMs within the 
sphere to have hits. 

The level 5 filter consists of containment cuts. These cuts 
are necessary as most of the remaining background events 
are located at the edges of the detector. The first contain- 
ment cut is on the reconstructed vertex position. This cut 
requires that the vertex of the event be more than 50 m in- 
side the top and bottom of the detector, and inside the outer 
ring of strings (string containment). The second contain- 
ment cut is on the charge received by each DOM. This cut 
requires that the DOM with the largest charge be located 
on an inner string (DOM charge containment). The effect 
of these containment cuts is shown in Figure 2. After these 
cuts are applied, background and signal simulation are fed 
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Figure 2: IceCube-40 detector in xy coordinates, black dots are the string positions. Each figure shows the reconstructed 
vertex for simulated E~ 2 signal, a) Before containment, b) after string containment, c) after DOM charge containment. 



into a machine learning algorithm. Multivariate analysis 
(TMVA) [12] is used, which assigns a boosted decision 
tree (BDT) response score to each event. This BDT re- 
sponse has a strong separation power. The variables used 
for machine learning are: 

• Z vertex position Reconstructed event depth. 

• Zenith track direction Reconstructed zenith angle. 

• Log likelihood Likelihood event arises from a muon. 

• Linefit velocity Particle speed to create hit pattern. 

• Eigenvalue ratio Topology of event. 

• Fill ratio Ratio of DOMs hit in sphere around vertex. 

• Time split Difference in time of two halves of event. 

• Containment Event distance from detector centre. 

The level 6 filter is the final level of cuts. There are two cuts 
at this level: the BDT response from multivariate analysis, 
and the reconstructed energy. These cut values are opti- 
mised using the Feldman-Cousins [13] method. The BDT 
response cut is 0.2, and the energy cut is 25 TeV. 

3.2 Filter Levels of High Energy Analysis 

This analysis was designed to quickly reach a result, and 
thus accepted a higher energy threshold [11]. It follows 
the mid energy analysis to level 3, then a BDT was trained 
using ten variables connected to shape, fit quality and po- 
sition of the event. Each variable has a correlation of less 
than 30% to reconstructed energy, so BDT score and en- 
ergy can be used as quasi-independent variables to suppress 
background in the final level. These cuts were optimised 
for sensitivity [14], based only on simulation. The back- 
ground is very sensitive to cosmic ray composition, ice and 
detector properties. In addition, limited statistics for rare 
events made the background prediction fragile. The atmo- 
spheric muon prediction is 0.72±0.28(stat) ±J || (sys) 
events, with a large systematic error that is based on a de- 
tailed comparison of simulation with data. The system- 
atic uncertainty will be reduced with future simulation. 
The signal prediction for an all flavor E~ 2 flux of 10 -7 
GeVsr^s^cm- 2 is 7.93±0.13(stat) ± 1.47(sys) events 
(assuming a flavor ratio of 1:1:1). 



4 Results 

4.1 High Energy Analysis 

In the high energy analysis four events were found. After 
careful inspection, all events look similar to background 
from atmospheric muons. Due to the systematic error this 
result is compatible with a null hypothesis. A limit at 90% 
confidence level on the astrophysical neutrino flux was set 
using the method of Rolke et. al. [15]: 



< 9.5 x 10- 8 GeVsr _1 s 



''cm 2 . 



(1) 



The energy range containing 90% of the signal is 89 TeV 
to 21 PeV. A comparison with model predictions and other 
analyses limits is shown in Figure 3. 
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Figure 3: Limits and flux predictions for an all flavor dif- 
fuse flux. Black lines indicate limits, this high energy anal- 
ysis is the solid line. The current best limit is given by 
the IceCube-40 diffuse search using muon neutrinos [16]. 
Theoretical models are thin grey lines. 



4.2 Mid Energy Analysis 

In the mid energy analysis a total of 14 events were ob- 
served. Of these, four events contained early hits, with 
timing from before the cascade. This could indicate a back- 
ground event due to an atmospheric muon interaction, a 
muon neutrino interaction, or muon production within the 



S. HlCKFORD etal. SEARCH FOR ASTROPHYSICAL NEUTRINO-INDUCED CASCADES 



cascade. The background prediction is under investigation. 
The remaining events contain no evidence of early hits, and 
after visual inspection appear to be high quality cascade 
candidate events, including three that have reconstructed 
energies above 100 TeV. 
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Figure 4: Distributions with final cut values shown by the 
vertical lines, a) BDT response, b) Energy spectrum. 




Figure 5: Cascade event in IceCube-40. Small dots are 
position of DOMs, circles show light received, where size 
indicates number of photons. 



events, or muons from muon neutrino interaction or within 
the cascade. Using atmospheric neutrino models Bartol and 
Sarcevic, 1.8 conventional events and 2.1 prompt events 
are predicted from atmospheric neutrinos. The events are 
currently being investigated by performing more extensive 
background simulations. 
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Abstract: We present the expected baseline sensitivity of the IceCube detector to cosmogenic neutrinos produced 
through the GZK process. Data from the partially completed IceCube detector have previously been searched for such 
highly energetic (> 10 6 GeV) neutrinos. With the completion of the detector in December 2010 and the full operation 
having started in May 201 1, IceCube's sensitivity to these neutrinos is significantly improved from previous studies. We 
calculate the expected sensitivity in the search of cosmogenic (GZK) neutrinos using a Monte Carlo simulation of the 
completed IceCube detector and the selection criteria developed in the previous analysis. The sensitivity for a diffuse 
flux of cosmic neutrinos with an E~ 2 spectrum in the central 90% energy range 300 TeV to 2 EeV is expected to be at a 
level of E 2 <f> Ve+ „^ +VT < 1.3 x 10 -8 GeVcm -2 sec _1 sr _1 with one year of operation. The corresponding differential 
sensitivity is also presented. 
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1 Cosmogenic neutrinos with IceCube 

Cosmogenic neutrinos are produced in the interactions of 
the highest energy cosmic -rays with the cosmic-microwave 
background (CMB) photons (the GZK process [1,2]) and 
subsequent charged pion decays [3]. These cosmogenic 
(GZK) neutrinos are one of the most promising messengers 
from the high energy, distant universe beyond PeV ener- 
gies. They may provide us with direct evidence of the high- 
est energy cosmic ray sources unlike the other messengers, 
such as gamma-rays and cosmic -rays, which experience in- 
teractions with the CMB and/or galactic and extra-galactic 
magnetic fields. 

IceCube is a cubic kilometer scale deep underground 
Cherenkov neutrino detector at the South Pole. The Ice- 
Cube detector construction was completed in December 
2010. The IceCube array [4] comprises 5160 optical sen- 
sors on 86 cables, called strings, over a 1 km 3 fiducial vol- 
ume of ice at a depth of 1450 m - 2450 m. In 2008-2009, 
40 out of 86 cables were deployed and taking data with 
an approximate fiducial volume of 0.5 km 3 . Results from 
the cosmogenic neutrino search with the half-completed 
configuration of IceCube [5] generated the best published 
limit to date on the neutrino fluxes above 1 PeV and up to 
10 EeV. 

In this proceeding, we present the expected sensitivity of 
the completed IceCube detector to cosmogenic neutrino 



fluxes calculated using a Monte Carlo simulation. The sig- 
nal discrimination methods are based on the selection crite- 
ria utilized in the 2008-2009 data analysis with the partially 
instrumented 40-string detector [5]. 

In the energy region above 1 PeV, the primary variable used 
to discriminate signal from background is the energy of 

zM 




Figure 1: Simulated single muon events in IceCube. Left 
panel shows a 100 TeV muon track representing the con- 
ventional event while the right panel indicates a 9 EeV 
muon EHE event. Circles denote optical sensors with more 
than one photo-electron signal recorded. The size of the 
circles represents the number of photo-electrons. Axes are 
distances in meters from the center of the IceCube detector 
array. 
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the particles. This is because the conventional atmospheric 
neutrino and muon background spectra are proportional to 
£-3.7 or s t ee p er5 w hile signal spectra follow E _1 <~ E~ 2 
in the energy region considered. Since the amount of en- 
ergy deposited in the form of Cherenkov photons by the 
neutrino-induced charged particles in the detector is highly 
correlated with their energy, the extremely-high energy 
neutrino signal stands out against the atmospheric muon 
and neutrino background because of the much higher light 
deposition. The total number of photo-electrons (NPE) 
recorded in an event is used as the main distinctive feature 
to separate signal from background. Figure 1 illustrates the 
difference in the energy deposition in the IceCube detector 
from a background-like 100 TeV muon and a signal-like 
9 EeV muon. 

2 Event selection 

The primary background in this analysis is muon bundles 
made up of large numbers of muons produced by high en- 
ergy cosmic -ray interactions in the atmosphere dominat- 
ing the downward-going directions. Because of the high 
multiplicity number, these events also leave a large amount 
of Cherenkov photons in IceCube. This background was 
simulated with the CORSIKA air-shower simulation pack- 
age version 6.720 [6] with the SIBYLL 2.1 [7] hadronic 
interaction model. Cosmic-ray interactions assuming pure 
proton and iron primary compositions in the energy region 
between 10 5 and 10 11 GeV were simulated. EHE neutrino 
signal events with energies between 10 5 and 10 11 GeV 
from several flux models were simulated using the JULleT 
package [12]. The cosmogenic neutrino induced tracks are 
most likely to have a near horizontal slightly downward- 
going geometry with falling distributions towards both ver- 
tically upward-going and downward-going directions due 
to the neutrino absorption in the Earth. 

The simulated high energy events are divided into the shal- 
low and deep event samples to take the difference in the 
optical properties of ice into account. The "depth" of the 
event is defined by the vertical position of the brightest 
photo-electron signal. The final background discrimina- 
tion is performed using different sets of variables for the 
shallow and deep events as described in Ref. [5]. Figure 
2 shows the event distributions in the planes of cos 9 vs 
NPE for the shallow events and AtLN-E vs NPE for the 
deep events. Here 9 is the reconstructed zenith angle of the 
event and AtLN-E is the time interval between the earli- 
est (E) and the brightest (LN, largest NPE) photo-electron 
signal in the event. A clear separation between the sig- 
nal and background can be observed. Reference [5] fur- 
ther describes the variables and compares the experimental 
and simulated event distributions. The straight lines and 
the quarter-elliptical shape show the applied NPE thresh- 
old value as a function of cos 9 and AtLN-E, respectively. 
The boundaries are set such that the background expecta- 
tion from cosmic-rays of an assumed pure iron primary is 
0. 1 events per year. For a pure proton case the background 



events are estimated to be at least a factor of 5 reduced 
from the current estimate of the background event numbers. 
This selection enhances the discovery potential of IceCube, 
which with a signal-to-background ratio of around 10 be- 
comes quite robust against large unknown systematics un- 
certainties in the background estimate. 

3 IceCube sensitivity beyond a PeV 

The quasi-differential model-independent sensitivity of the 
IceCube detector at 90% CL per energy decade for neu- 
trino fluxes above 10 15 eV (1 PeV) is shown in Fig. 3 
assuming full standard neutrino mixing. The correspond- 
ing sensitivity for a diffuse flux of cosmic neutrinos with 
an E~ 2 spectrum in the central 90% energy range from 
300 TeV to 2 EeV is calculated to be E 2 ^ c+ ^ + „ T < 
1.3 x 10 -8 GeVcm -2 sec~ 1 sr _1 with one year of observa- 
tion. The improvement of the sensitivity from the analysis 
of the data taken with the half completed IceCube [5] is 
approximately a factor of two. 

Table 1 gives the event rates for several model fluxes of 
cosmogenic neutrinos assuming cosmic -rays to be protons 
only. We expect 0.8 to 1.7 cosmogenic neutrino events 
per year, assuming moderate to strong cosmological source 
evolution models, while 0.11 background events are ex- 
pected in the same time period. 

The corresponding neutrino effective area is shown in 
Fig. 4. The neutrino effective area represents the surface 




log NPE log NPE 



Figure 2: Event number distributions of the shallow (upper 
panels) and deep (lower panels) event samples in 365 days 
are shown for signal (left panels) and background (right 
panels) simulations. The signal distributions are from the 
cosmogenic neutrino model in Ref. [8] adding all three fla- 
vors of neutrinos. The background distributions are from 
CORSIKA-SIBYLL with iron primaries. The lines in each 
panel show the final selection criteria. 
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Figure 3: All flavor neutrino flux quasi-differential sensitivities of the IceCube detector after one year (filled squares) 
and five years (filled stars) of operations. Several model predictions (assuming primary protons) are shown for compar- 
ison: YT [8], Kalashev et al. (strong evolution) [9], ESS (Qa = 0.7) [10], Ahlers et al. (maximal), Ahlers et al. (the 
best fit, incorporating the Fermi-LAT bound) [11]. The gray dashed horizontal line indicates the Waxman-Bahcall flux 
bound with cosmological evolution [17]. Model fluxes are summed over all neutrino flavors, assuming standard neutrino 
oscillations. The model independent differential upper limits by other experiments are also shown for Auger (PAO) [13], 
RICE [14], ANITA [15], the previous IceCube results (2007-2008, IC22) [16], and (2008-2009, IC40) [5]. Limits from 
other experiments are converted to the all flavor limit assuming standard neutrino oscillation and a 90% CL per energy 
decade quasi-differential limit when necessary. 



area of an equivalent detector if it were 100% efficient. For 
v„ and v T , the areas exceed 10 3 m 2 at 10 9 GeV which is 
the main energy range in the IceCube cosmogenic neutrino 
search [5]. The present analysis is sensitive to all three 
neutrino flavors. Similarly the effective area near the Ice- 
Cube detector can be defined as the area within which the 
neutrino induced muons and taus, or neutrinos are 100% 
detectable. They are shown in the right panel in Fig. 4. Ice- 
Cube acts as a detector with effectively 50% larger volume 
than its actual size for neutrino-induced muons at 10 9 GeV. 
The effective area near the detector for the neutrinos inter- 
acting near or inside the detector (direct neutrino channel) 
is more than two orders of magnitude smaller than those for 
muons and taus. However, the direct neutrino interactions 
form an important contribution because the neutrino flux at 
IceCube depths is two orders of magnitude larger than the 
flux of the secondary charged leptons. 

4 Discussions 

IceCube may be the first experiment to probe the cosmo- 
logical evolution of the cosmic-ray sources [19]. 
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Figure 4: The left panel shows 47r solid angle averaged 
neutrino effective areas for each neutrino flavor. The 
dashed line corresponds to v e +D e . The solid line is v^+D^, 
and the dotted line is v T + P T . All assume equal flux of neu- 
trinos and anti-neutrinos. The right panel shows the effec- 
tive detection area near the detector for secondary muons, 
taus and primary all flavor neutrinos. 

While models of astronomical neutrinos include uncer- 
tainties in the photon field the cosmic rays interact with 
prior to escape from sources, CMB induced cosmogenic 
neutrino models are not affected by this uncertainty. In- 
stead, these neutrino fluxes are highly dependent on the 
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Kalashev et al. [9] 


5 
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1.65 ±0.01 


Ahlers et al. dip transit at 10 19 eV (best) [11] 
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2 
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0.80 ±0.01 


Ahlers et al. dip transit at 10 19 eV (max.) [11] 


4.4 
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10 21 eV 


1.69 ±0.01 



Table 1: Expected numbers of events by IceCube in 365 days from several cosmogenic neutrino models assuming the 
cosmic-ray primaries to be protons. The spectral indices 7, cutoff energies E max at sources as well as cosmological 
evolution indices to and extensions in redshift Z max for the cosmic-ray sources are also listed for reference. The corre- 
sponding expected number of background events in one year is 0.11 ± 0.01. Errors are statistical only. 



cosmological distributions of the cosmic -ray sources, the 
cosmic-ray energy spectra in the sources, and the cosmic- 
ray composition. However, the cosmogenic neutrino event 
rates expected in IceCube are relatively stable under differ- 
ent assumptions on the spectra injected at the cosmic-ray 
sources, such as the maximum proton energy at sources and 
spectral indexes. Reference [9] shows that the neutrino flux 
significantly increases with decreasing primary spectral in- 
dex from 2 to 1 only in the energy region above 10 18 GeV. 
It is also shown in Ref. [10] that when the spectral index 
increases to 3, the neutrino spectrum is shifted to a slightly 
lower energy region. For both cases, only a slight increase 
in the event rates is expected since the target neutrino en- 
ergies of IceCube are below 10 185 eV [5]. Similarly the 
dependence on E max is weak because it mainly affects the 
flux shape above 10 18 GeV [9]. Therefore IceCube is sen- 
sitive to the redshift evolution of the co-moving density of 
sources, often parametrized with an exponential index to 
as (1 + z) m , to include the redshift dependence. 

Studies of the baseline capability of the IceCube detector 
for cosmogenic neutrinos in the energy region above 1 PeV 
have shown that IceCube is able to detect cosmogenic neu- 
trinos or constrain fluxes with moderate parameters, m = 
3.5 - 4. 

In this paper we only considered the case of pure proton 
cosmic-ray primaries for GZK neutrino production. A de- 
tector with an order of magnitude larger fiducial volume 
than IceCube, such as the next generation radio Cherenkov 
detector arrays, ARA [20] and ARLANNA [21], is re- 
quired for measurements of neutrino fluxes induced by 
iron-dominated cosmic -rays. 

This is a performance study based on the method already 
realized and applied to the actual experimental data anal- 
ysis [5]. New techniques are being developed for fu- 
ture analyses [22, 23]. For example, there is a tech- 
nique to experimentally identify atmospheric muon bun- 
dle events by looking for a coincidence signature in the 
IceCube optical modules and the IceTop air shower sur- 
face array [22]. Differences in the energy loss profiles and 
waveform shapes between the atmospheric muon bundles 
and neutrino-induced single muon events are also being ex- 
plored. These new techniques should improve sensitivities 
for future cosmogenic neutrino analyses with IceCube. 
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Abstract: The IceCube neutrino telescope was constructed to search for high energy neutrinos of cosmic origin. At 
the highest energies, neutrinos associated with the interaction of the most energetic cosmic rays with cosmic microwave 
background photons (GZK effect) are considered a guaranteed signal, with expected event rates of up to a few events 
per year in a cubic kilometer detector. Searches for GZK neutrinos have been performed using data taken with the 
intermediate construction stages of the now complete IceCube detector. We present the results of finished and on-going 
analyses, with a focus on the search using data taken between spring 2009 and spring 2010, when the IceCube detector 
was roughly 70% complete. 
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1 Introduction 

The detection of extremely-high energy (EHE) neutrinos 
with energies in excess of 10 7 GeV may shed light on the 
yet unknown origin of the highest energy cosmic rays. 
The direct observation of EHE charged cosmic ray parti- 
cles is limited by their inevitable energy loss in the cos- 
mic microwave background through photo-pion produc- 
tion, known as the Greisen-Zatsepin-Kuzmin (GZK) ef- 
fect [1]. The trajectories of the charged cosmic ray par- 
ticles with diminished energies will have been random- 
ized in cosmic magnetic fields upon arrival at the Earth. 
Neutrinos from the decays of the secondary charged pions, 
tt^~ — > ^Vfj, — > e ± v e v li v ll , will travel in straight lines and 
unattenuated over cosmological distances and carry infor- 
mation about the sources of EHE cosmic rays. 

The IceCube neutrino observatory consists of a cubic kilo- 
meter sized Cherenkov detector embedded in the 2800 m 
thick glacial ice cap at the South Pole and an overlying 
square kilometer surface air-shower array. The in-ice de- 
tector consists of 5160 light sensitive digital optical mod- 
ules (DOMs) deployed at depths between 1450 and 2450 m 
on 86 vertical cables ("strings"). Each DOM is equipped 
with a 25 cm photo-multiplier tube (PMT) along with two 
waveform digitizers and supporting data acquisition, cali- 
bration, and control hardware [2, 3]. Interactions of high 
energy neutrinos with the surrounding matter are detected 
via their Cherenkov emissions in the highly transparent Po- 
lar ice [4]. With its large detection volume, the in-ice de- 



tector is the first neutrino telescope with a realistic chance 
to detect the small flux of EHE neutrinos associated with 
the GZK effect. 

During IceCube's construction phase, which started in 
2005, data taken with the partially instrumented in-ice neu- 
trino telescope have been searched for signatures of EHE 
neutrinos [5, 6]. The analysis of data taken during the years 
2008/2009, when 40 of the 86 strings of the in-ice detector 
were deployed, has led to the currently most stringent lim- 
its on fluxes of EHE neutrinos with energies between 1 PeV 
and 10 EeV (Figure 1). In these proceedings, we report on 
a search for EHE neutrinos in data taken with the 59-string 
detector between spring 2009 and spring 2010. 

2 Method 

The vast majority of the events recorded by IceCube are 
due to down-going atmospheric muons that are created 
by interactions of high energy cosmic rays in the atmo- 
sphere, and which are sufficiently energetic to penetrate 
the ice overburden and deposit Cherenkov light in the de- 
tector. Against this background, an EHE neutrino interac- 
tion inside or in the vicinity of the detector would stand 
out with a much higher Cherenkov light deposition. Fig- 
ure 2 shows the expected light deposition in terms of the 
number of recorded photo-electrons (NPE) and its corre- 
lation with the zenith angle (cos 6) of the primary parti- 
cle tracks for simulated GZK neutrino induced events and 
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Figure 1 : Present flux limits on EHE neutrinos compared 
to two all flavor GZK neutrino flux predictions, YT ((m, 
Z max ) = (4,4) [7] and ESS (fi A = 0.7) [8], and the 
Waxman-Bahcall bound [9]. Limits from IceCube 40- 
string GZK neutrino search (IC-40 GZK, differential limit 
and E~ 2 spectrum integrated limit) [6] are shown as thick 
black lines, the limit from an UHE neutrino search (IC-40 
UHE) [10] is shown as grey horizontal line. Less stringent 
limits were set by Auger (PAO) [11], ANITA [12], and the 
IceCube 22-string detector (IC-22) [5]. 

simulated atmospheric background events. For both signal 
and background events, the light yield is strongly correlated 
with cos 0. For down-going atmospheric muon events, the 
expected NPE rapidly decreases with increasing inclina- 
tion of primary particle track, because with increasing slant 
depth the muons will lose more energy before reaching the 
detector. Below the horizon (cos < 0), low energy at- 
mospheric neutrinos arriving from the opposite hemisphere 
are the only expected background. The detection probabil- 
ity of GZK neutrinos is highest for directions close to the 
horizon, because of the relatively short neutrino interaction 
length at EHE energies. 

Exploiting the correlation between the event light yield 
and the track direction, the data selection criteria to sep- 
arate signal from background events in IceCube's EHE 
neutrino searches have routinely been designed as two- 
dimensional boundaries in the NPE-cos 0-plane [5, 6]. 
Simple geometric fit methods, whose performance proved 
robust against systematic uncertainties in the detector re- 
sponse, were used to infer the track directions. In this anal- 
ysis, we follow the same strategy, but we use a different 
fitting algorithm than previous analyses to determine the 
track direction. 

Following a blind analysis procedure, the selection criteria 
are optimized on simulated signal and background events. 
A subset of 10% of the experimental data, evenly dis- 
tributed throughout the data taking period, is used to val- 
idate the detector simulation. After the selection criteria 
are developed, the data selection is applied to the blinded 
90% of the data, which for the 2009/2010 data-taking pe- 



riod roughly comprises 330 days of detector livetime. The 
10% subset is discarded, in order to avoid statistical bias. 



3 Monte Carlo simulations 

The dominant background at the final data selection lev- 
els is high multiplicity muon bundles induced predomi- 
nantly by heavy cosmic ray primaries with PeV to EeV 
energies. This background was simulated with the COR- 
SIKA air-shower simulation [13] using the SIBYLL 2.1 [14] 
hadronic interaction model. Two primary types, proton 
and iron, with energies between 10 4 and 10 11 GeV were 
simulated. The primaries were sampled from a power-law 
energy spectrum following dN/dE oc E~ 2 , in oder to 
over-sample the high energy end of the cosmic ray spec- 
trum, which is most important to this analysis. Proton and 
iron components are then re- weighted to broken power-law 
spectra, whose combination approximates the all particle 
spectrum at PeV energies and above [15]. 

Signal events induced by EHE neutrinos in the energy 
range between 10 5 and 10 11 GeV were simulated with 
the JULleT package [16]. The charged secondary particles 
were sampled from an energy spectrum oc E ~ 1 . The events 
can be re-weighted to various GZK neutrino flux predic- 
tions. In these proceedings, we use the predictions from 
references [7] and [8] (c. / YT and ESS in Figure 1). The 
quoted event rates correspond to the sum of all three neu- 
trino flavors, v e , v^, and v T . 

4 Event selection 

The first data selection for the various IceCube analyses is 
performed on-line at the South Pole, before data are sent to 
the northern hemisphere. For this analysis, the on-line filter 
required a minimum of 630 photo-electrons to be recorded 
in an event. 

Following the analysis strategy that was developed for the 
40-string detector [6], further data reduction is achieved by 
requiring at least 200 DOMs to have registered light within 
a time window of [— 4.4/is, 6.4/zs] around the largest local 
light deposition in the detector. The latter is defined as the 
time at which 10% of the largest PMT pulse was captured. 
Further, we require the total number of photo-electrons 
recorded in this time window to be larger than 3200. With 
these requirements, the atmospheric background is reduced 
by two orders of magnitude, while 75% of the signal is re- 
tained (Table 1). While the previous EHE search used the 
linefit algorithm [17, 5] to reconstruct the track directions, 
we use a different algorithm in this analysis, the dipole- 
fit [17]. The dipolefit assigns a dipole moment, M, to the 
light pattern recorded in each event: 



M = 



1 



2 



-Ndom 

E • 

; £dx>m +1 
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Figure 2: Expected light yield (logarithm of the number of detected photoelectrons, log 10 (NPE)), for simulated signal 
and background events and a subset of the experimental data set (left) and the correlation with the cosine of the zenith 
angle of the primary particle tracks for simulated atmospheric background events (middle) and GZK neutrinos (right). 
The background prediction includes cosmic ray induced muons, assuming a cosmic ray spectrum according to [15], and 
atmospheric neutrinos according to [18] with a contribution from prompt neutrinos according to [19]. The GZK neutrino 
spectrum is simulated according to [7] (c. / YT in Figure 1). 



cut Experimental Background GZK YT [7] GZK ESS [8] 

NPE > 630 6.20 x 10 7 (6.86 ± 0.18) x 10 7 2.35 ± 0.01 1.81 ±0.01 

NPE > 3200 and N UO m > 200 6.65 x 10 5 (7.68 ± 0.12) x 10 5 1.80 ±0.01 1.38 ±0.01 

log 10 (jVP£)-0.5-:P>4 336 365 ±6 1.39 ±0.01 1.07 ±0.01 



Table 1: Expected event rates at various selection levels for 330 days of detector livetime. The signal rates correspond to 
the GZK neutrino models according to [7] (YT) and [8] (ESS). Errors are statistical only. 



where fj is the vector of spatial coordinates of the DOM 
that recorded the ith light signal in time, and A^dom is the 
total number of DOMs fired in the event. The magnitude of 
the dipole moment \M\ takes values between and 1, and 
provides a measure for the directionality of the light flow in 
the event: large values of \M\ indicate a track-like signal, 
while small values indicate a rather spherical light pattern. 
EHE neutrino interactions typically yield small dipole mo- 
ments. Cascades induced by v e and v T interactions natu- 
rally generate spherical light patterns, and the light pattern 
from induced EHE muons has a broad radial distribu- 
tion. Low energy atmospheric muon events on the other 
hand, typically have dipole moments close to 1 . 

Figure 3 shows the correlation between the magnitude 
and the direction (cos m) of the dipole moment for at- 
mospheric background events and GZK neutrino signal. 
Background events cluster in the region with large dipole 
moments and down-going directions. A combination of 
both magnitude and direction of the dipole moment, V = 
cos 9m + 2 • |M|, is used as a measure of the similarity 
of an event to a down-going track. Compared to a cut on 
the reconstructed direction only, a larger fraction of the ex- 
tremely bright EHE neutrino signal events is preserved. 

Background events with low values of V are predomi- 
nantly induced by low energy neutrinos and atmospheric 
muons that pass outside the instrumented volume and de- 
posit only very little Cherenkov light. Figure 4 shows 
the correlation of V with the light yield NPE. Signal 
and background events are well separated in the V-NPE- 



plane. A two-dimensional cut in this plane defined by 
(\og 10 (NPE) — 0.5 • T>) > 4 reduces the atmospheric 
background by two more orders of magnitude, while the 
expected GZK neutrino signal still exceeds one event in 
330 days (Table 1). 

Additional selection criteria to separate the GZK signal 
events from the remaining background are being investi- 
gated. A realistic chance to detect a GZK neutrino signal 
requires the further selection criteria to keep the signal ex- 
pectation above 1 event, while suppressing the background 
to a level of 0(0.1) expected events per year. The GZK 
neutrino search with the 40-string detector [6] achieved 
a signal expectation of 0.5 events for GZK models pre- 
sented here, above an expected background of roughly 0. 1 
events at the final selection level. With the larger 59-string 
detector, improved event selection criteria, and better un- 
derstanding of the detector response, the required signal 
to background ratio to either detect GZK neutrinos or to 
constrain the here considered flux predictions seems within 
reach. 

5 Conclusions 

The detection of GZK neutrinos with IceCube seems tan- 
talizingly close. The analysis of data taken with an inter- 
mediate construction stage of the detector, in which half 
the in-ice detector was deployed, allowed to place the most 
stringent limits on EHE neutrinos to-date. Data taken with 
later construction stages are presently being analyzed. Al- 



IceCube EHE neutrino search 




Figure 3: Observables of the dipolefit. The magnitude of the dipole moment | M\ versus the reconstructed zenith angle 
(cos 6m) for simulated background events (left) and simulated signal events (middle). The linear combination V = 
cos 9m + 2 • \M\ (right), is a measure the similarity of the event to a down-going track. 




log (NPE) log (NPE) log (NPE) -0.5© 



Figure 4: Correlation of the combined dipole observable V = cos0 M + 2 • \M\ with NPE for simulated background 
events (left) and simulated signal events (middle). The linear combination of T> and NPE that is used as a cut parameter 
is shown in the right panel. The separating boundary defined by log W (NPE) — 0.5 • V > 4 is shown as black lines in 
the two dimensional plots. 



ready the next construction stage, which roughly 70% of 
the detector components deployed, may reach the sensitiv- 
ity to probe current models of GZK neutrino fluxes. 
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Abstract: The detection of cosmogenic (GZK) neutrinos with IceCube requires the ability to discriminate very rare and 
energetic signal events from an abundant background of cosmic ray induced muons. High energy cosmic ray air showers 
produce high numbers of muons densely packed around the shower core trajectory. These bundles of muons emit large 
amounts of Cherenkov light in the ice that constitutes the detection volume. We present several techniques to improve 
background rejection while keeping a large fraction of the GZK neutrino signal. The differences in the light distributions 
around a neutrino-induced muon track and a muon bundle are exploited. The photon hit-time pattern in the detector 
differs slightly for the two event types and is used for identification of muon bundles. The surface array, IceTop, is used 
to tag the background with high efficiency but limited zenith range. The efficiency of this method was studied using data 
from the partially completed detector. 
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1 Introduction 

Ultra High Energy Cosmic Rays (UHECRs) with en- 
ergies above 10 11 GeV have been observed by several 
experiments[l, 2]. The origin of UHECRs remains un- 
known, though there may be indications of a correlation 
of incoming directions with the close-by extra-galactic 
source distribution [3]. The elucidation of the origin has 
been longed for from the first detection. UHECRs interact 
with cosmic microwave background photons and necessar- 
ily generate neutrinos with energies in excess of 10 7 GeV 
through secondary pion decays (GZK effect). Therefore, 
the detection of such Extremely High Energy (EHE) neu- 
trinos can shed light on the UHECR origin. 

The IceCube detector [4], completed in Dec. 2010, instru- 
ments a huge volume of 1 km 3 ultra transparent glacial ice 
and is suitable to search for rare EHE neutrino events. 

The GZK neutrino flux prediction depends on the cosmo- 
logical source evolution, the source injection spectra and 
the cosmic ray composition [5]. The expected GZK neu- 
trino event rate in IceCube is about one event per year 
[5, 6]. 

The main background for EHE neutrino search comes 
from muon bundles induced by cosmic ray interactions in 
the atmosphere. While bundles are much more abundant 



compared to the expected neutrino signal, their flux de- 
creases steeply with increasing energy. Therefore signal, 
expected to have a harder energy spectrum, may emerge 
from the background above a certain critical energy. In ad- 
dition, muon multiplicity in bundles increases with the pri- 
mary cosmic ray energy, which leads to more pronounced 
background-event signatures and to increasing rejection ef- 
ficiency. Another difference between neutrinos and muon 
bundles is their arrival direction. While the muon bundle 
rate decreases with increasing zenith angle, near horizon- 
tal directions are favored for GZK events because of the 
increase of the neutrino cross section at high energies. 

The energy and arrival direction information has been used 
in several EHE neutrino searches [7, 8] producing the best 
upper limit for EHE neutrinos in the relevant energy range 
around 10 9 GeV. 

In this paper we present methods which are being de- 
veloped to achieve higher signal efficiency and high- 
multiplicity muon bundle background rejection using the 
characteristic differences between the two. 

2 Muon Bundle Rejection Techniques 

In Extensive Air Showers (EAS), more than thousands of 
muons can be generated and reach the IceCube detector 
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Figure 1: Top: Average number of recorded pulses as 
a function of the time residual for events with 10,000 
<NPE< 30,000. Bottom: negative time residual likelihood 
ratio distribution for the same events. Dashed red: Vu sig- 
nal simulation, dotted blue: background simulation, black: 
experimental data. 

depth, under an ice overburden of about 1500 m. Most 
of these muons are concentrated in a dense core, but some 
may have relatively high transverse momenta p t and are 
therefore separated from the core of the bundle at the depth 
of IceCube by a distance oc pt/E^. Multiple scattering and 
deflection due to the Earth's magnetic field can increase 
the separation for near-horizontal events [9]. Observing the 
separation of single muons within the bundle core is not 
possible in IceCube as the photon scattering length in ice 
is too short and the detector's string spacing is too large for 
this purpose. However, differences in the light distribution 
around the bundle core compared to that around a single 
muon can be used to distinguish the two event classes. 

2.1 Early Photon Hit Times 

The application of a single muon hypothesis track recon- 
struction [10] to a muon bundle event gives the location and 
direction of the bundle axis. For each detected Digital Op- 
tical Module (DOM) pulse and a given reconstructed track 
we define the time residual t ICS as the difference between 
the measured pulse time and the expected arrival time of 
an unscattered Cherenkov photon from the single track hy- 
pothesis. In the muon bundle case, the light generated by 
outlying muons may result in pulses with negative t res val- 
ues, indicating photon arrival times inconsistent with the 
single track hypothesis. We exploit the density of nega- 
tive t ICS pulse distribution by means of a likelihood analysis 
where signal and muon bundle background hypotheses are 
compared. The distribution of number of pulses with nega- 
tive t ICS values for simulated signal and background, com- 
pared to experimental data events for the IceCube 40 string 
configuration is shown in Fig. 1 (top). The observable NPE 
refers to the total Number of Photo-Electrons collected in 



Figure 2: Top: Distribution of detected light perpendicu- 
lar to a simulated neutrino-induced single muon track (red 
markers) and the fitted /i function (black) for an event with 
100,000 NPE. Bottom: light distribution likelihood param- 
eter ratio distribution for events with 10,000 < NPE < 
30,000. 

an event by the IceCube DOMs. The resulting likelihood 
parameter distributions are shown in Fig. 1 (bottom). 

2.2 Perpendicular Light Distribution 

The amount of detected light at perpendicular distances 
from a single muon track is a function of the muon energy, 
ice properties and the detector noise level due to DOM 
electronics. Its parameterization is described in [11] and 
is given by 

fi(d,e i E)=a(9,E)u- d / do +b noisc (1) 

with d: DOM-track distance, 9: string axis-track opening 
angle, E: energy of the track and do — 1 m. Parameter 
a (in units of NPE) represents the light normalization and 
is dependent on energy of the track in IceCube, the dimen- 
sionless parameter ui describes the shape of the falling light 
curve and 6 no i S e (in units of NPE) gives the expected noise 
level of the DOMs. Parameters a and oj are both dependent 
on ice properties which vary with depth [12]. However, 
it is difficult to resolve the dependency as for each event 
light is emitted and detected at different ice depths and the 
dependency is averaged in the fitted parameter values. 

Examples for the detected light distribution around a sin- 
gle muon neutrino event and the fitted /i(d, 9, E) function 
are given in Fig. 2 (top). The radial spread of the muons 
inside the bundle can be up to ~ 50 m, so the perpendic- 
ular light distribution at small distances is flatter around 
a bundle compared to a single muon track. The amount of 
detected light at larger distances is higher for a single muon 
track compared to a bundle. This could be because muons 
in the bundle range out and do not reach the clearest ice at 
the bottom of the detector and because of large stochastic 
energy losses in the single muon case. These differences 
appear mostly in the value of the fitted to parameter. For 
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Figure 3: Combined A ln(L) cut signal efficiency vs. back- 
ground efficiency. For a signal cut efficiency of 90% the 
background rejection efficiency ranges from 90%-98% for 
neutrino energies of 10 6 5 — 10 11 GeV. The curves refer to 
samples of events with increasing amount of light detected 
with the IceCube 40 string detector. 



Figure 4: effective areas for a simple analysis using a 
combined A ln(L) cut on MC simulation of IceCube 40 
string configuration. Full red: IceCube EHE events fil- 
ter level, dotted purple: MRF cut calculated on filter level 
event sample, dashed blue: MRF cut calculated after ap- 
plying a combined A ln(L) cut. 



signal and background events in a GZK neutrino search the 
obtained 2-dimensional distributions of fitted uj versus a 
values occupy different areas of the phase space. A likeli- 
hood parameter comparing signal and muon bundle back- 
ground hypotheses is constructed and its distributions for 
the IceCube 40 string configuration simulated and experi- 
mental data events are shown in Fig. 2 (bottom). 

2.3 Results 

The bundle rejection observables described in sections 2.1 
and 2.2 were combined in a single likelihood parameter. 
Background vs signal efficiency is shown in Fig. 3 as a 
function of a cut on the combined likelihood ratio A ln(L) 
defined as ln(L s i g ) — ln(Lbg)- In order to assess the 
strength of the combined A ln(L) observable, a cut was set 
at a fixed A ln(L) value which gives signal passing rates of 
88% — 97% depending on NPE range. A signal selection 
NPE threshold was then calculated using the MRF tech- 
nique [14] on simulated and experimental data that passed 
both the EHE event filter (NPE > 630) and the combined 
A \n(L) cut. The resulting effective areas are shown in 
Fig. 4. 

3 IceTop Veto on Cosmic Ray Showers 

An EAS event in the IceCube detector may be preceded 
by hits recorded in the surface detector IceTop. Therefore 
another promising technique to discriminate muon bundles 
from EHE neutrinos is to use IceTop to veto muon bundles. 
IceTop uses the same DOMs as IceCube to detect the elec- 
tromagnetic and muonic part of EAS. For the IceTop veto 
the electromagnetic part plays only a minor role. Inclined 
EAS are mostly tagged by detection of high p t muons far 
away from the shower core. 



An efficient IceTop veto against high energy EAS will im- 
prove the signal efficiency in a GZK neutrino search. The 
three main parameters determining the EAS veto efficiency 
are: 

Primary Energy and Composition: The higher the pri- 
mary energy of the EAS, the higher the probability to see a 
signal in IceTop, as the number of secondary particles and 
the lateral extension increases. A higher veto probability 
for heavier primary particles is expected due to more sec- 
ondaries. 

Distance to IceTop: The shorter the distance of the shower 
core from IceTop, the higher the probability to detect the 
event by IceTop. The distance of the shower core can be 
up to several km, depending on the geometrical hit posi- 
tion in IceCube. This parameter is closely related to the 
inclination. 

Inclination: With increasing inclination the air shower 
propagates through more atmosphere where the electro- 
magnetic shower component gets attenuated more than the 
hadronic component. Thus for near-horizontal showers we 
expect IceTop to detect mainly muons. 

Single tank hits in IceTop are used to establish the IceTop 
veto [13]. These IceTop hits have to be within the time 
window of several /is of a high energy event that triggered 
IceCube. In order to find hits in coincidence with the air 
shower front, we reconstruct the shower front in time and 
space. Here the center of gravity of the IceCube event and 
the direction from the track reconstruction are used. A pla- 
nar shower front is a good approximation to find coinci- 
dent hits and can be corrected by a parameterization of the 
shower front curvature. Fig. 5 shows the principal idea of 
the IceTop veto. 

The distribution in Fig. 5 is used to fix the size of the veto 
time window to 400 ns covering the coincidence peak. For 
comparison, we take a background time window of the 
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Figure 5: Sketch illustrating the veto principle. The right 
plot shows the IceTop hit distribution over time relative to 
the shower front. The sharp peak at T = in the plot is 
caused mainly by coincident IceTop hits. The second peak 
after the shower front has passed is caused by after-pulses 
of the photomultipliers. 

same length before the shower front reaches IceCube so 
that only hits from uncorrelated cosmic ray showers are in- 
cluded. The upper plot of Fig. 6 shows the time window 
selection strategy. The lower plot shows the number of 
events as a function of hits for both the random and real co- 
incidence time windows. The ratio of the number of events 
in the real coincidence time window over the number of 
events in the random coincidence time window for a given 
number of hits gives the probability of the event to have 
correlated hits. Due to the low statistics of the available 
test data (10% of IceCube data taken in 2010 with 79 string 
configuration), no detailed systematic analysis of the veto 
efficiency is possible yet. The random coincident back- 
ground is following Poisson statistics and is independent 
of the coincident hits. We assume that our data contains 
exclusively cosmic ray induced events and neglect all other 
background contributions. Estimates of the veto efficiency 
have to take into account the full dependences on energy, 
inclination, and distance to IceTop. All events with an NPE 
value exceeding 10 5 in the IceCube 2010 test data fulfill the 
99.99% coincidence probability (5 hits in the coincidence 
time window). We conclude that the high veto efficiency 
at very high energies makes the IceTop veto a powerful in- 
strument for GZK and diffuse neutrino searches. 



4 Possible Future Improvements 

Other variables are under investigation with the aim of fur- 
ther improving the sensitivity to EHE neutrinos. The longi- 
tudinal distribution of the amount of detected photons has 
been found to differ between a single muon derived from 
a neutrino and muon bundles [15]. The observed distribu- 
tion along muon bundles is rather smooth while the one 
around single muons fluctuates much more and exhibits 
DOMs that detect very low PE values, which do not exist in 
the bundle case. It is important to separate EHE neutrinos 
passing far away from the detector center from the abun- 
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Figure 6: Top: Illustration of the chosen coincidence time 
windows for events with NPE > 30,000. The background 
time window contains hits due to background from uncor- 
related cosmic ray showers. Bottom: The number of events 
in the background time window falls steeply as a function 
of IceTop hits. The distribution for events in the coinci- 
dence window is significantly flatter, indicating a prelim- 
inary veto efficiency higher than 85% for NPE > 30,000. 



dant lower energy muon bundles passing well within the 
fiducial detector volume since both of them yield similar 
NPE. Time-over-threshold of the recorded charge for each 
DOM was found to be a good proxy for the distance of light 
source from the DOM. This information will be used in a 
future EHE neutrino search as well as utilizing the IceTop 
veto power. 
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